Phase diffusion in graphene-based Josephson junctions 
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We report on graphene-based Josephson junctions with contacts made from lead. The high 
transition temperature of this superconductor allows us to observe the supercurrent branch at tem- 
peratures up to ~ 2 K, at which point we can detect a small, but non-zero, resistance. We attribute 
this resistance to the phase diffusion mechanism, which has not been yet identified in graphene. By 
measuring the resistance as a function of temperature and gate voltage, we can further characterize 
the nature of electromagnetic environment and dissipation in our samples. 

PACS numbers: 74.45,+c, 74.50.+r, 73.23.-b, 72.80.Vp 



Josephson junctions with a normal metal region 
sandwiched between two superconductors are known 
as superconductor-normal-superconductor (SNS) struc- 
tures. Over the years, the normal region has been 
made from no n- metallic nanostructures, including het- 
erostructures, nanotubes, quantum wires, quantum dots 
(H, and, most recently, graphene Usually, these 

superconductor-graphene-superconductor (SGS) junc- 
tions employ aluminum as the superconducting metal, 
separated from graphene by another metal layer (often ti- 
tanium) intended to create a good contact. In this paper, 
we succeed in making palladium- lead (Pd/Pb) contacts 
to graphene. Here, Pd is known to form low-resistance 
contacts to graphene 0, 0] , while Pb has the advantage 
of a relatively large critical temperature (7.2 K). As a 
result, the SGS junctions demonstrate an enhanced zero- 
bias conductance up to temperatures of the order of 5 
K, and at temperatures below ~ 2 K a clearly visible 
supercurrent branch appears in the I — V curves. 

In all of our samples, a small, but non-zero voltage 
is observed below the switching current. We attribute 
this feature to the phase diffusion mechanism [8]. The 
phase diffusion in underdamped junctions is enabled by 
the junction's environment, which provides dissipation 
at high frequencies [9[ . Observation of this regime in our 
SGS junctions is facilitated by the high critical temper- 
ature of Pb. We first study the phase diffusion resis- 
tance as a function of temperature, which allows us to 
extract the activation energy associated with the phase 
slips. Next, the phase diffusion is measured at differ- 
ent gate voltages, resulting in a consistent picture of 
the junction's environment and dissipation at high fre- 
quencies. This series of measurements allows us both 
to establish the phase diffusion regime in underdamped 
SGS junctions, and to analyze their behavior in terms of 
well-established models. Finally, we demonstrate an effi- 
cient way of controlling the junction by passing a current 
through one of the electrodes within the same structure: 
the locally created magnetic field modulates the critical 
current. Several periods of oscillations are visible, indi- 
cating the spatial uniformity of the junction. 

Graphene was prepared by a version of the con- 



ventional exfoliation recipe [To| from natural graphite 
stamped on RCA-cleaned Si/Si02 substrates. The sam- 
ples were verified by Raman spectroscopy to be single 
atomic layer thick with low defect density The 
electrodes were patterned by standard e-beam lithogra- 
phy and thermal evaporation. We first deposited ~ 2 
nm of Pd, which formed highly transparent contacts to 
graphene [6|, |7j, followed by ~ 70 nm of Pb. Care was 
taken not to heat the samples above ~ 90° C, and to store 
them in vacuum in order to minimize oxidation of Pb. 

The inset of Figure 1(a) shows a scanning electron 
micrograph of a typical device. A layer of graphene is 
visible as a gray triangular-shaped shadow in the cen- 
ter of the image, contacted by two long metallic elec- 
trodes. A known current (AC + DC) is driven through 
the graphene between two probes on one side of the sam- 
ple, and voltage is measured between two probes on the 
other side. We present the results measured on three 
different samples. Sample A has a gap of d = 100 nm 
between the leads; the graphene region is L = 1.5/im 
long. In samples B and C, the leads meander across 
graphene for a much longer total distance of L ~ 15 jam 
and ~ 20/im, respectively (see schematics in Figure 3). 
The gap between the leads is designed to be d = 500 nm 
(B) and 400 nm (C). 

Figure 1(a) demonstrates the simultaneously measured 
DC voltage V and differential resistance dV/dl vs. ap- 
plied current / in sample A. (The inset shows a different 
sample of a similar design.) From the dV/dl curves, it 
is clear that a pronounced effect of superconductivity is 
observed at temperatures as high as ~ 5 K, which is 
comparable to the transition temperature of the leads 
(verified to be « 7 K). At the two lowest temperatures, 
the I — V curves show a region of vanishing small V; the 
junction abruptly switches to a normal state when the 
current exceeds a certain value (the switching current, 
Is)- On the reverse current sweep, voltage drops close to 
zero at the retrapping current (Ir). Figure 1(b) plots Is 
and Ir at the two lowest temperatures vs. the normal 
resistance of the sample, controlled by V ga t e - 

Observation of the hysteresis in the I — V curves in- 
dicates that the junction is underdamped [8]. Indeed, 
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FIG. 1: a) Top inset: scanning electron micrograph of a typ- 
ical sample and the measurement schematic. Two Pd/Pb 
contacts are made to graphene (gray triangular shade). A 
fixed DC current / with a small AC modulation (tens of nA) 
is driven through graphene between contacts on one end of 
the sample, and the voltage drop V is measured between two 
contacts on the other end. Bottom inset: schematic showing 
the sample layout and defining dimensions. Main panel: DC 
voltage V and differential resistance dV/dl vs. bias current I 
measured at several temperatures on sample A (different from 
the sample shown in the inset). V ga te = +40 V is applied to 
enhance the conductance of graphene. Each curve is mea- 
sured while sweeping the current from negative to positive, 
resulting in hysteresis at the lower temperatures, at which a 
difference appears between the switching and retrapping cur- 
rents. (The spikes in dV/dl at the switching and retrapping 
currents are naturally truncated in the measurement.) b) 
Switching and retrapping currents (Is and Ir) as a function 
of Rn, which is controlled by the gate voltage. The normal 
resistance is extracted from the I — V curves as dV/dl at a 
current of exceeding the switching current; thus defined 
Rn virtually does not depend on temperature. Stars: critical 
current extracted as I a = cEa/Fl from the activation energy 
Ea of phase diffusion (see Figure 2 for more details). Inset: 
maps of dl/dV vs. I and V ga te at 5 different temperatures. 



the estimated quality factors of our junctions are of the 
order of one (see also the discussion of Figure 3). Here, 
we take into account the presence of the degenerately 
doped Si substrate, which provides the dominant contri- 
bution to the capacitance between the superconducting 
leads (tens of fF). An alternative explanation of hystere- 
sis in a SNS junction could be overheating [l2[. In our 
case, two samples (A and C) have very similar switch- 
ing and retrapping currents. Their normal resistances, 
which control the heat generation just before the retrap- 
ping, are different only by a factor of ~ 2. However, the 
dimensions of graphene regions, which control the heat 
dissipation, are vastly different: the areas differ by ~ 50, 
and the contact lengths differ by ~ 15. Therefore, con- 
ventional underdamping, rather than overheating, seems 
more likely in our case. 

In all our samples, a finite voltage on a /iV scale ap- 
pears on the superconducting branch of the I — V curve. 
This behavior is illustrated in Figure 2a, showing the 
I—V curves measured in sample B at three different tem- 
peratures, including 1.4 K, at which the I — V curve is 
hysteretic. The appearance of a finite voltage is explained 
by the "phase diffusion" mechanism, where a point rep- 
resenting the phase slowly descends the tilted washboard 
potential [8[, getting trapped at successive local minima 
following each phase slip. The existence of the phase dif- 
fusion regime in an underdamped junction indicates an 
efficient high-frequency dissipation due to the junction 
environment [9]. Experimentally, we find that the mea- 
sured values of the switching current are reproducible 
upon successive sweeps, again supporting the phase diffu- 
sion mechanism as opposed to premature switching by a 
single phase slip While not yet reported in graphene, 
the phase diffusion regime has been recently analyzed in 
a conceptually similar case of an underdamped junction 
based on a multiwall carbon nanotube 
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The presence of phase diffusion allows us to investigate 
the rate of phase slips, proportional to the sample resis- 
tance, and its dependence on temperature. Theoretically, 
the zero-current differential resistance due to the phase 
diffusion should depend on temperature as [q|, [l3 - [l7| 



Ro(T) (x T- x ex V (-2Ejlk B T). 



(1) 



Here, the Josephson energy Ej = HI^ /2e, and 1^ is the 
true critical current of the junction. Figure 2(b) shows 
the product TRq(T) plotted as a function of the inverse 
temperature for samples A and B. Both samples clearly 
show activation behavior; the extracted activation energy 
turns out to be close to twice the Josephson energy, as 
estimated from the switching current. This is illustrated 
in Figure lb, where for the ease of comparison we con- 
vert the activation energy Ea to current as I a = cEa/Fl 
(stars), which is indeed close to Is- 

Let us now analyze the dependence of the phase diffu- 
sion resistance Rq on Ej, controlled by V ga te- Let us 
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FIG. 2: a) I — V characteristics of sample B at several tem- 
peratures and Vgate = 0. Finite voltage could be noticed 
below the switching current at the lowest temperature, b) 
The product of the temperature times the differential resis- 
tance, TRo, as a function of inverse temperature 1/T, mea- 
sured on sample A (filled symbols), and on sample B (empty 
symbols). In sample B, V ga te = 0, while in sample A, several 
values of V ga te are taken, resulting in several sets of symbols. 
Evidently, in all sets, TRo demonstrates activation behavior, 
with an activation energy of Ea ~ 10 K. This energy is con- 
verted to critical current according to I a = eEA/fr, shown 
by stars in Figure 1(b). (The differential resistance at small 
current, Ro(T), becomes too small to measure at low Rn, 
so the analysis is limited to the high- R n range.) c) Sym- 
bols: the prefactor to the exponential, R' (see text), vs. Ej 
extracted from the same data as in panel (b). Lines are a 
linear fit, assuming R f ~ ZqEj /JzbT, which corresponds to a 
junction underdamped at DC but overdamped at the plasma 
frequency. For comparison, the dashed lines illustrate the ex- 
pression Rq oc VEj, which clearly does not fit the data well, 
d) The product of the switching current and the normal re- 
sistance IsRn vs. inverse resistance 1/Rn- 



define the prefactor to the exponential in eq. (1) as 
R' = R e 2Ej / kBT . Theoretically, this prefactor varies 
depending on whether the Josephson junction is over- 
damped or underdamped. For an overdamped junction, 
R ~ REj/k B T 



15 



16j, where R is the shunting re- 
sistance, i.e. ~ R N . In case of an underdamped junc- 
tion, Rq ~ -^hujp /ksT [9j, so that R' depends on Ej 
and the junction capacitance C through the plasma fre- 
quency ujp oc y^Ej/C. Finally, if the junction is under- 
damped at DC, but overdamped at the plasma frequency, 
Rq scales as ex ZqEj /kpT where Zq is the real part of 
the environmental impedance at high frequency [17j . 

Since C and Zq do not change with the gate voltage, 
while Rn and Ej do, we may distinguish between the 
different cases. In Figure 2(c), we plot R = R§e 2Ej l kBT 
vs. Ej (taken as Ea/2) for three temperatures, 1.3, 1.7, 
and 3.0 K. It is clear that the scaling of R' is consistent 
with oc Ej and is not consistent with either R f oc ^/Ej 
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FIG. 3: (a) Map of the voltage drop V measured vs. bias cur- 
rent / and the current II, which flows along one of the leads 
parallel to the interface with graphene and induces magnetic 
field Bl- The dark regions along the horizontal axis corre- 
spond to the supercurrent branch. The current is swept from 
the negative to the positive direction, resulting in the visible 
hysteresis between the retrapping (negative I) and switching 
(positive /) currents in the central lobe, (b) Extracted Is and 

Ir VS. I L . T = 1.3 K, Vgate = 40 V. 



or Rq oc RnEj (not shown) [18]. This observation al- 
lows us to identify the junction as underdamped at DC, 
with plasma frequency oscillations damped by the envi- 
ronment; the environmental impedance is found to be 
Zq « 200 — 250 tt. The overall agreement convinces 
us that the macroscopic behavior of the junction is ade- 
quately described by Ref. 

Using Ref. LL7| we estimate that at T = 1.3 K, Is is 
close to I(p (exceeds 70% for the whole range shown in 
Figure 1(b)). Therefore, we can use Is in place of 1^ 
and plot IsRn vs. 1 /Rn in Figure 2(d). The trend in 
the graph resembles that of 1^ Rn vs. the Thouless en- 
ergy, Eth-> as expected in the SNS junctions [l9j]- Indeed, 
Eth should be inversely proportional to the resistivity of 
graphene. At the location of the "knee" in the curve, Eth 
is estimated to be of the order of A, indicating the transi- 
tion between the ballistic and diffusive SNS regimes. We 
do not attempt a more careful comparison of these pre- 
liminary data with theory, since extracting Eth from Rn 
would require the exact knowledge of the contact resis- 
tance and the density of states in the sample. Also, the 
superconducting gap is likely suppressed at the interface, 
which would complicate analysis. 

Let us now discuss the effects of magnetic field on the 
junctions. To generate the field, we passed a large (mA 
range) DC current II along one of the Pb leads, parallel 
to the interface with graphene. For these measurements, 
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we picked relatively large pieces of graphene, and made 
the junction's length L tens of /am by meandering the 
leads across the sample surface (samples B and C). The 
resultant large area between the leads allowed us to pass 
several flux quanta through graphene, before II drove 
the Pb lead normal [2o| . 

Figure 3a shows the DC voltage drop across graphene 
V, mapped as a function of the current / flowing through 
graphene, and the current II generating the magnetic 
field. Several regions of vanishing voltage are visible 
along the horizontal axis. From the extent of these re- 
gions, one can extract the switching and retrapping cur- 
rents, Is and Ir, vs. Il (Figure 3(b)). The resulting 
modulations are close to the expected Fraunhofer pat- 
tern / ex sin(7r/L/^L,o)/^L [s|, where Il,o corresponds 
to passing one flux quantum through the junction. Ob- 
servation of several oscillations (about 5 at both positive 
and negative iz,, not shown) indicates a uniform junction. 
We also found that similar modulation are induced by an 
externally applied magnetic field, as reported previously 
in other S-graphene-S samples [2|45j]. 

The difference between Is and Ir, which exists in the 
center of the pattern in Figure 3a, disappears at higher 
magnetic field, i.e. for lower Is- For example, it is 
not seen in the side lobes at all (Figure 3b). This im- 
plies that the quality factor Q(4° } ) = (2el$ ] '/hf^R 
is close to 1 at the central lobe. Assuming that other pa- 
rameters of the junction, except for Iq\ do not depend 
on magnetic field, we may fit Ir as f(Q)I^\ where is 
f(Q) is a universal function, approximated at Q ~ 1 as 
f (Q) « 1.273 - 0.311Q - 0.030Q 2 + 0.013Q 3 [21]. We 
can further replace 1^ with the measured Is (see e.g. 
(22| ) - indeed, based on Ref. the two currents are es- 
timated to be very close for Is > 0.4/iA at T = 1.3 K, as 
we have already discussed for sample A. The fit shown in 
Figure 3(b) is achieved by taking Q = 1.4 at the center 
of the pattern as the only fitting parameter. As expected 
from the theory [2l|, the difference between Ir and Is 
disappears at Q « 0.85. 

In conclusion, we describe a simple method of making 
S-graphene-S Josephson junctions that operate at tem- 
peratures of up to several Kelvin. All of our samples 
demonstrate phase diffusion - a small, but finite differ- 
ential resistance at zero current, with an activation en- 
ergy close to twice the Josephson energy. We analyze 
the phase diffusion in some detail, and find it in good 
agreement with the established theory for a junction un- 
derdamped at low frequencies and overdamped at the 
plasma frequency. We also demonstrate the efficient con- 
trol of the critical current flowing through graphene by 
running a current through one of the leads within the 
same structure. The observed modulation pattern shows 
several periods, indicating the spatial uniformity of the 
junction. 
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